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We study the Friedmann-Robertson- Walker model with dynamical dark energy modelled in terms 
of the equation of state p x — w x (a(z))p x in which the coefficient w x is parameterized by the scale 
factor a or redshift z. We use methods of qualitative analysis of differential equation to investigate 
the space of all admissible solutions for all initial conditions on the two-dimensional phase plane. 
We show advantages of representing this dynamics as a motion of a particle in the one-dimensional 
potential V(a). One of the features of this reduction is the possibility of investigating how typical 
big rip singularities are in the future evolution of the model. The properties of potential function V 
can serve as a tool for qualitative classification of all evolution paths. Some important features like 
resolution of the acceleration problem can be simply visualized as domains on the phase plane. Then 
one is able to see how large is the class of solutions (labelled by the inset of the initial conditions) 
leading to the desired property. 
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^ ■ T. INTRODUCTION 

^sO \ While the cosmic acceleration becomes one of the most exciting discoveries of the cosmology |2j the nature of the 
dark energywhich is still unknown 0, 0] . Unfortunately, there are very weak constraints on its form of an equation 
of state [2, l|. The most popular candidate for dark energy is the cosmological constant A, which can be treated as 
some kind of perfect fluid satisfying the equation of state p = —p, p = A. However, it remains to explain why value of 
A, obtained from type la supernovae (SNIa) observations, is so small in comparison with the value of A interpreted 
as vacuum energy (i.e., the Planck mass scale). Recent investigations in this field concentrate mainly on the dark 
energy, modelled as (a) quintessence scalar field 6] or (b) based on barotropic equation of state ■ In this paper we 
q-( assume the general form of the equation of state in which the equation of state factor w of total matter depends on 
I ■ redshift z through the scale factor a, so w(a{z)) = p/p, w < —1/3 in some nonempty interval of evolution. 

It is characteristic the existence of large variety of different dark energy models [||. We unify all these models 
assuming that dark energy obeys the general form of the equation of state p = w(a(z))p. It is interesting to find some 
general properties of an evolutional path of such models. For this aim dynamical system methods seem to be a natural 
method because it offers the possibility of investigating whole space of solutions starting from all admissible initial 
conditions. The philosophy of qualitative investigation of differential equations shifts interest from founding their 
exact solutions toward investigating their sensitivity (or fragility) with respect to small changes of initial conditions 
or model parameters. 

It is demonstrated that dynamics of the FRW model with dark energy can be represented in the form of two- 
dimensional dynamical system, x — P(x,y), y — Q{x,y) where P,Q G C°°. The phase portraits of this system are 
organized by critical points (xq, yo) : P(xq, yo) — Q{xq, yo) = or limit cycles (non-point-like attractors). Following the 
Hartman-Grobman theorem in neighborhood of non-hyperbolic critical points (ReXt ^ 0, where are eigenvalues 
of linearization matrix x = Ax, A is the Jacobi matrix) the nonlinear matrix is equivalent to its linearization. The full 
knowledge of dynamical behavior requires analyzing its behavior at infinity. To achieve this it is useful to transform 
trajectories from a R 2 phase plane into a Poincare sphere Then infinitely distant points of the plane are mapped 
into a sphere's equator S . The type of critical points are conserved under this mapping but critical points present in 
infinity can appear at the equator. Hence an orthogonal projection of any hemisphere onto the tangent plane gives a 
compactificd phase portrait. 

The main subject of this paper is presentation how dynamics of the FRW model with dynamical dark energy can 
be reduced to the form of a two-dimensional dynamical system (section [HJ). The general properties of such a system 
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are investigated in section ITTT1 We find that a global structure of the compactified phase space RP 2 depends on the 
strong energy condition only. In section HVI we present advantages of using a complementary description of dynamics 
as a Hamiltonian flow. In this picture the evolution of the model is represented as a motion of unit mass in the 
one-dimensional potential V(a). It is interesting that a shape of the diagram of V(a) contain all informations which 
are needed to determine all critical points at the finite domain as well as at infinity and their character. The detailed 
discussion of this issue and the compactified phase portraits contain section^ We find only three generic cases of the 
phase portraits. All FRW models with a postulated form of dark energy belong to the one from distinguished class 
of phase portraits modulo homeomorphism preserving orientation of the phase curves (which establish equivalence of 
the phase portraits). Therefore it is possible to classify all FRW dark energy models on the phase portraits without 
some details about the specific form of the equation of state. 



II. THE FRW MODELS WITH DYNAMICAL DARK ENERGY AS A DYNAMICAL SYSTEMS 

Let us consider the FRW model with source in the form of noninteracting dust matter and dark energy which the 
equation of state is parameterized by p x = w x (a)p x , w x < —1/3. Then from the conservation condition 

p = -3-(p+p), (1) 
a 

where a dot denotes the differentiation with respect to time t, we obtain for each component relations 

Prn = Pmfi a 3 > Pm(o> = 1) = Pm,Q] 



Px = p x ,oa 3 exp i 



Px (a=l)= Px . a . (2) 

Hence we can combine both effects and define the total pressure as 

p = + w x p x = w(a)p, (3) 

where the equation of state factor for total matter is 

i \ w x (a) . 

W W = n —th r~\ > ( 4 ) 

1 + S^expfs/ ^Mda 



where fi^o are density parameters for both matter and dark energy, a is expressed in the units of its present value ao. 

In our further analysis we assume the equation of state in form ©. Note that the coefficient of state w x (a) can be 
always extracted from the integral differential equation 

w x (a) = iu(a)jl + ^y^exp ~~~ da ) j i ( 5 ) 



equivalent to the differential equation 



d y- s^.ijHW, y _™M 



da a w(a) 

which the solution, for the initial condition w(l)/w x (l) = ^ix,o/(^x,o + ^m.o)j is in the form 

_ w x (a) n m /_ f a w x (a) 



w(a) 



= l + ^exp(3 r^da-l). (6) 
Six,o v Ji a J 



Due to © it is possible to determine w x (a) provided that w(a) is given. 
The FRW dynamics is governed by equations 

k a 2 
p = 3^ + 3^, 

'a, a 2 k 

P=- 2 2 "I ' 



(7) 
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where p = p c g and p are effective energy density and pressure respectiveiy; p c ff = Pm + Px in the model under 
consideration, k - the curvature index k = 0, ±1. System {7J) together with condition (JIJ is closed if the form of the 
equation of state © is postulated. 

Then the basic dynamical equations reduce to 

a = y, 

y = -\{p{a) + ip(a)}a. (8) 
6 

Note that expression 

£ + V(a)«-£, (9) 

plays the role of first integral (JHJ if 

V(a) = -Pcffy. (10) 

Equations JHJ constitute a two-dimensional autonomous dynamical system which solutions can be visualized in the 
phase plane (a,y). Then 10 represents the equation for an algebraic curve on which lies solutions of the dynamical 
system. Among all solutions there exist singular solutions represented by critical points in the phase plane. At the 
finite domain of the phase plane they always represent static universes (a = 0) 

j/o = and p(a ) + 3p(a ) = 0. (11) 

They are located on a- axis, and exist if only p < —p/3 for effective pressure and energy density. There are two 
possibilities: p(ao) = or w(ao) — —1/3 for appearing such a points. If w(a) is a monotonic function of a then there 
is one solution ao : w(ao) = —1/3. Note that for all critical points the function V(ciq) = — k/2. Because it is always 
the non-positive function if only p c s > the critical points are admissible in the region occupied by trajectories of 
the closed model (fc = +1). 

Following the first integral © the whole phase space (a, y) is divided on two disjoint domains by the trajectory of 
the flat model. In these regions we find closed models in the region y 2 /2 + V(a) < while open models are situated 
in the region y 2 /2 + V(a) > of the phase space. The boundary of the strong energy condition p + 3p = defines 
the line in the phase space which separates regions occupied by non-accelerating and accelerating phases of evolution 
of the models. 



III. THE GENERAL PROPERTIES OF THE FRW MODELS WITH DYNAMICAL DARK ENERGY 

The assumed form of the equation of state includes different dark energy models. For example the Chaplygin gas 
can be treated as a special case of some prescribed form of the equation of state ■ In this case we have 

V = = w(a)p, (12) 

P 

where w(a) = — Aa 6 / (Aa 6 + B) and A, B are positive constants. From (|12fl one can seen that w varies in the interval 
[—1,0]. It is a monotonic (decreasing) function of a. 

If we consider the equation of state for a non-interacting mixture of dust (p = 0) and radiation matter then we 
obtain the equation of state in the form 

1 1 Pr 1 1 

p = Q+oPr=o — ; P=o o P> ( 13 ) 

3 3 p r + p m 3 "m,0 

o a 

"r,0 



which can be reproduced by w — w(a) 
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Freese and Lewis [t| have recently considered models (called Cardassian) which are alternative to dark energy models. 
Following their approach the source of acceleration is in some modification of FRW dynamics instead of dark energy 
conditions. One can show that the Cardassian models are equivalent to dark energy models for two components of 
matter. The basic equations are 



PcS = (Pm + Pr) + B(p m + p r ) n , 

3H 2 = p eS , 

where H — (In a)' is the Hubble function and p r is energy density of radiation matter: k 
model, n is the model parameter and 



(15) 



is the assumption of the 



PeS = W(a)(p m + Pr) + w(d)(p m + p r ) 7 



(16) 



where w(a) = n(w + 1) — 1 and w(a) is given by formula (|14fl . Therefore, assumed for © of equation of state seems 
to be sufficiently general to unify all cosmological models which offers explanation of SNIa data (note that it is also 
valid of the Linder and Copeland model) . 

Because we assume the presence of dust matter we have xo ^ but in general for the fluid violating the weak 
energy condition (p + p){x) > it is possible to have the critical point xq = yo = 0. Then det A = l/6(p + 3^)1(0,0) 
and A 2 + det A = 0, therefore eigenvalues of the linearization matrix are real of positive signs which correspond to 
saddles points. Note that in this case there are no static critical points at p + 3p = in the finite domain. 

The phase space of the model is organized by critical points and phase curves determined from the first integral 
@ . The critical points at the finite domain of the phase space are situated on the a-axis as intersection points of the 
a-axis with the boundary curve of the strong energy condition p + 3p > 0. The character (type) of critical points is 
determined from eigenvalues of the linearization matrix A: 
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■~{xp(x)(l+3w(x))) 



(17) 



J {xo,yo) 



where x = a is a dimensionless scale factor expressed in terms of its present value ag. 

The eigenproblem det [A — Al] = then reads A 2 — XTiA + det A = 0. Consequently the sign of the determinant A 
determines the type of the critical points, i.e., whether A is real or complex. It is consequence of the fact that 



TrA = and 
A 2 + det A = 0. 



detA=^(xp(x)(l + 3w(x))) 



(18) 



Because at the critical points we have (p + 3p)(xo) =0we have 



, , 9 

det A = xo 7t — 
ox 



(p + 3p)(x) 



(19) 



The critical points of the system can be saddle points if det A\ Xo < or centers in the opposite case if det A\ Xo > 0. In 
the first case eigenvalues are real of opposite signs while in the second case they are purely imaginary and conjugated. 
The determinant of the linearization matrix can by rewritten to the form 



det A = 



dp(x) 



dx 



(l + 3w(x)) 



+ 3p(x a ) — 
ax 



(20) 



Therefore if we consider critical points appearing at w(xq) = —1/3 then det A = 3p(xo)dw /dx\ Xo and the type of 
critical points depends on sign of dw/dx at the critical point xq (we assume p{x) > 0). If w(x) is a decreasing function 
of the scale factor then the critical points are saddle types. 

If we consider the second type of critical points which is the static critical point located at x = xo such that p(xo) = 
then first term in H20JI decides whether the critical points are saddles or centers. Therefore if xq : w(xq) < —1/3 
(x : w(x ) > —1/3) then det A < (det A > 0) if dp/dx\ Xo > (like for phantom fields violating the weak energy 
condition p + p > 0). If xq is located in the decelerating region Xq ■ w(xq) > —1/3 then det A < (det A > 0) if 
dp/dx\ xo < (like for matter satisfying a weak energy condition). If det A > or det A < 0, centers or saddles are 
admissible respectively. Note that both energy density of matter and dark energy are positive this case is excluded at 
finite domains. Finally one can conclude that if we assume positivity of energy of each component then the existence 



5 



as well as character of the critical points are determined by energy conditions in the neighborhood of the transition 
moment corresponding the inflection of the diagram a(t). From the physical point of view it is moment of hanging 
decelerating phase into accelerating one. 

The full knowledge of dynamical behavior requires an analysis of the system at infinity. One can perform such an 
analysis in a simple way by introducing the projective coordinates on the phase plane. There are two maps which 
cover a circle S* 1 at infinity 

1 y 

(1) (z,u): z——, u — —. z = 0, — oo < u < +oo, 

1 X x (21) 

(2) (v,w): v = — , w = — , v = 0, — oo < w < +oo. 

V V 

By adding a circle S 1 at infinity to R 2 we obtain a compact space - projective plane HP 2 . In the (z, u) coordinates 
the system under consideration has the following form 

dz 

"77 = - z «> 

7 (22) 
du 1 2 

-77 = -g(p + 3 P )- U , 

with the first integral in the form u 2 + 2z 2 V(l/z) = kz 2 . There are two types of critical points situated on the circle 
at infinity (z — 0,— oo < v < +oo) 



2 = 0, u = ±J~(p + 3p) . (23) 

Of course such points are admissible if the strong energy condition at the critical point is violated (p + 3p)o < and 
the value of (p + 3p)o is finite (we can also find a critical point at a finite domain which was detected previously 
uq = (H = oo) and (p + 3p)o — 0). From the form of the first integral of 12211 we obtain that critical points 123|) 
are intersection points of the trajectory of the flat model u 2 + 2z 2 V(l/z) = with the circle at infinity. They are 
representing stationary solutions Hq = ±y/—l/6(p + 3p)o, i.e., the expanding and contracting deSitter solutions. For 
existence of such points it is required that z 2 V(l/z) goes to a constant value as z — > 0. For example if we consider the 
ACDM model with the cosmological term A then p(z = 0) = A, Hq = ±^/l/3. In the special case if z 2 V(l/z) — > 
as z — > we obtain a static solution. 

For the analysis of the type of critical points we consider the linearization matrix 

-Uq 

l(0,u ) 



In this case TrA = — 3uq which means that if uo > (expanding deS+) the deSitter space is a global attractor, while 
if uq < deS_ is representing a global repellor. If uq = all eigenvalues Ai,2 = l/2(— 3u ± |tio|) are degenerated, 
non- hyperbolic critical point. 

Analogical investigations of the critical points at infinity can be performed in the map (v,w): v — l/y,w = x/y. 
The dynamical system (JSJ in these coordinates is of the form 

dv vw , „ . , W . 

dt & v (25) 

dw w , „ . , w . 

= 1+ p + 3p - , 
dt 6 v 

with first integral l/(2u 2 ) + V(w/v) = -k/2. 

At circle at infinity v = the effects of curvature are negligible. Therefore the critical points at infinity are 
intersections points of the algebraic curve 1 + 2v 2 V(w/v) = with a circle at infinity v — {y — oo). If the function 
V is negative and for large x behaves asymptotically like V{x) cx x m , i.e., it is a homogeneous function of degree m, 
then for m — deS attractors (repellors) are admissible. If m > 2 then there are two critical points in the domain 
y > and their counterparts in the symmetric domain y < (with respect to an x-axis). They lie on the trajectory 
of the flat model and a circle S 1 . For both critical points the projective coordinates are the same Vq = 0, Wq = 0. 
If we consider V(x) in the neighborhood of x = in a power law approximation: V(x) oc x m , m < then the 
critical point x = 0, y = oo (v — 0) can be recovered. In the neighborhood of this critical point, representing an 
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initial singularity, dust matter effects dominate while the effects of dark energy (for which w x < —1/3) are negligible. 
Big rip singularities arise in the phantom dark energy cosmological model • At a big rip type of singularities the 
cosmological scale factor and the Hubble function achieve an infinite value in a finite interval of time. For both critical 
points at infinity y = ±00 (also H — 00) but for the big-rip singularity (in contrast to big bang when x = 0) we have 
x = 00 (or w = 0). If we put m = 2 + e, where e > then y oc x 1+£ approximate trajectories near big-rip singularity 
while yjx achieve infinite value at big-rip singularity as well as at the initial singularity the values of x are infinity 
and zero respectively, so they appear as a different point on the phase plane (x,y). 

Finally one can conclude that existence and type of critical points at the finite domain an at infinity is determined by 
the strong energy conditions (at finite domain) and the trajectory of the flat FRW model. The intersection trajectory 
of the flat model with a circle at infinity determines the location of the critical point at infinity. We assume a two- 
components model of matter filling the Universe. In the neighborhood of the initial singularity x = dust matter 
dominate the effects of dark energy. Therefore for the late time dark energy (w x < —1/3) should be feasted. Three 
different scenarios are possible if dynamical behavior for late time is admissible. It depends whether w x + 1^0 and 
consequently the sign of w x + 1 determines the type of a global attractor (a deSitter is w x = — 1 or a big-rip attractor 
if Wr < —1 or a model with asymptotic H = 0). From the form of first integral © one can observe symmetry of 
the phase space y — > —y. For the case of w x > — 1 there is no deSitter attractor but trajectories going toward the 
state ao : p(xo) — the static critical point at infinity. Note that while the localization of the static critical point and 
its character depends on solutions of the equation (p + 3p)(a) — and d/da| a= o(p + 3p) in the neighborhood of ao, 
the behavior at infinity (t — > 00) depends crucially on parameter m in asymptotic of V oc x m . 

Note that this ty pe of the future singularity may appear for the pressure (to = —00) and not only for the energy 
density p (see also 11]). Note that because our system is autonomous it has symmetry t — > t + b , b =const it is 
always possible that a singularity occurs for some finite value of the cosmological time. In other words any two phase 
portraits with rescaled time (through a smooth function) are the same. 



IV. ACCELERATING FRW UNIVERSE AS A PARTICLE IN A ONE-DIMENENSIONAL POTENTIAL 

In this section we adopt the Hamiltonian formalism to our dynamical problem. This gives at once the insight into 
the previously introduced function V(x). Moreover our problem stays similar to that of a particle moving in the 
one-dimensional potential energy V(x). To prove this let us rewrite Eq. (jSJ) to the form analogous to the Newtonian 
equation of motion 

dV\ 1 f x 

V(x) = - / (p + 3p)ada. (26) 



dx J 6 Jo 

If we assume that p(a) and p(a) satisfy the conservation condition Q the integration by part can be exactly 
performed. Then, after including Q rewritten to the equivalent form 

s<*>+4<*>-°- (27) 

which enables us to eliminate p from V given by H26[). we obtain, as a result of integration, 

V(a) = -p£, (28) 

where if a is expressed in the units ao, V(a) = V{x). 
Of course, system 1126^ has the first integral 

y + V(x) = ±Q ki0 , (29) 

where here a dot denotes differentiation with respect to rescaled t time r: t — * r: dt — \Hq\ _1 dr, flk,a — —k/a^ is 
the density parameter for curvature fluid, V(x) — — (l/2)f2 m oa;^ 1 — SJii 2 ; Q m ,o and fl x are density parameters for 
matter and dark energy, respectively. 

Therefore, from Eq. I|29|l we find interpretation of V(x) as a potential function for a particle-universe moving in 
the configuration space x: x > 0. Equation lf^|l can be simply reduced to the two-dimensional dynamical system 

x = y, 

dV 

V ~ "to' (30) 

1L + v(x) = \si k , Q . 
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In terms of the potential critical points are represented as extremal points: dV/dx\ Xo = 0. The linearization matrix 
and eigenvalues are 

A=( J^V I), (31) 

v dx 2 u y 



cP-V 

A 2 + Vd =0- (32) 
dx 2 lx=x o v ; 

From Eq.© we obtain that if d 2 V/dx 2 \ XQ < 0, i.e., a diagram of a potential function is upper convex, the 
eigenvalues are real of opposite signs which corresponds to saddle points. In the opposite case, d 2 V/dx 2 \ XQ > 0, 
eigenvalues are purely complex and conjugate which correspond to centers. The saddle points in contrast to centers are 
structurally unstable critical points. While the saddles appear at maxima of V(x), the centers correspond to minima. 
Therefore the diagram of the potential function with two maxima has to possess a minimum which corresponds to 
a center located between two saddles in the phase space. But the presence of the center in the phase space make 
the system structurally unstable. Hence the situation in which only one static critical point exists in finite domain is 
representing generic situation. The rest-exceptional cases are mere complicated and numerous and they interrupt the 
discussion of typical garden variety of cosmological dynamical systems. 

Early history of investigations in the theory of dynamical systems has been dominated by searching for generic 
properties. They define a class of phase portraits that are far simpler that arbitrary one [T^ | . The part of the Kupka- 
Smale theorem states the genericity of Gl for critical points. A dynamical system has the property of Gl if all of its 
critical points are hyperbolic (or elementary). The center is an example of a non- hyperbolic critical point. Therefore 
the corresponding dynamical systems are exceptional or non-generic, 

From eq. (|2(j[l we obtain that if V(x) is an increasing function of its argument then in the interval of a; a universe is 
decelerating, while in the interval of x in which V(x) is a decreasing function of x, a universe is accelerating. At the 
maximum the universe starts accelerating which corresponds to redshift zt : 1 + zt — Xq 1 . This type of dynamical 
behavior in the neighborhood of redshift transition is a prototype of evolutional scenario of the dark energy epoch. 



V. THREE GENERIC SCENARIOS OF COSMOLOGICAL DYNAMICS WITH DARK ENERGY 

The results of the qualitative analysis of the dynamical system are summarized and presented in Fig. While all 
phase portraits are topologically equivalent (modulo homeomorphism preserving an orientation of the phase curve) 
at the finite domain, some differences are in the location of the critical points at the circle at infinity. There are only 
three prototypes of dynamical behavior of the model. 

On the phase portraits we consider the plane (a, a) where the coordinate a is positive. All these solutions which 
correspond to the flat universe are located on the curve a 2 = — 2V(a). Trajectories going inside of two branches of 
this curve correspond to closed models while those moving outside correspond to open models. 

One can distinguish the vertical boundary line {a : dV/da — 0} which separates the solutions with acceleration 
from those without acceleration. The intersection of this boundary with the a-axis is at the saddle which is quite 
similar to Einstein's static universe. The cosmological singularities lie on the circle at infinity a 2 + a 2 = oo. They 
always appear as a intersection of a trajectory of the flat model and the circle at infinity. On the circle at infinity S* 1 
which compactifies the plane R 2 to the projective plane, the effects of curvature are negligible. 

Let us concentrate on the behavior of trajectories at the finite domain which is identical for all cases in Fig. ^CD We 
find only one static critical point ao located on the a-axis which corresponds to a maximum of the potential function. 
Trajectories moving in the region a < ao confined by separatrices approaching and escaping saddle (on the left from 
saddle point denoted as /) correspond to the closed universes expanding to the maximal size and then contracting to 
the final singularity. Note that they correspond to the closed universes which never undergo an acceleration phase 
during the whole evolution. There exists also two another types of evolutional paths for closed models. They are (we 
describe expanding trajectories) located between the upper branch of flat trajectory and the separatrices which lies 
in the region y > 0. They start from the standard singularity a = 0, a = oo then reach a point of maximal expansion 
d and continue expansion to the deS+ node (or state H = in the case in Fig. [3). 

Another type of evolution experienced by closed models lies in the region (IV) which boundaries form two sepa- 
ratrices, one coming from the saddle point to deS+ and other one coming out the singularity and approaching the 
saddle. The trajectories located in this region describe the evolution from the state a = oo, d = oo (d = in the 
exceptional case in Fig- EI) deS_, i.e., the contracting deSitter model; next they reach a point of minimal contraction 
d = and then they expand to deS+ model if V(a) oc a m and to = 2 or reach a state d = oo, a = oo if to > 2. In 
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other words the trajectories running in the region (/V) undergo a contraction from the unstable contracting deSitter 
node toward the stable deS + model - the expanding deSitter model only if m = 2 (i.e., the cosmological constant case 
H = const). In the other case there is no deS+ as a global attractor in the future (H = oo). It should be pointed 
out that the existence of deS+ attractor in the future is the unique property of the dark energy model for which 
V oc p x x 2 /6 and p x (oo) — const > 0. If to < 2 and V oc x m as x goes to infinity then trajectories start and land at 
the degenerated critical point (x = oo, x = 0) which is of course the non-generic critical point. It is interesting that 
the location as well as the type of the critical points at infinity can be determined from asymptotic behavior of the 
potential function as x — > oo (or x — > 0). 

Let us briefly comment the generic case when V(x) oc x 2+e , e > 0. Then in the future evolution of the model with 
e > instead of deS phase there is unwanted type of a singularity, so called a big-rip singularity. Because this critical 
point is a stable node big-rip singularities are generic. Note that this type of final behavior is a global attractor for all 
open, closed and flat models. All these models undergo the transition from the decelerating to accelerating phases. 
All these models also undergo the so-called loitering phase like for the Lemaitre model 01 . In the language of the 
potential function if the particle universe starts off with just enough kinetic energy to take it up to the top of the 
potential then value of kinetic energy close to this critical value gives rise to the loitering phase. Finally the existence 
of a loitering phase is a generic property of the models with dynamical dark energy. 

The phase of possible phase portraits which are generic can be classified by the construction of representative cases 
(Fig. n and |2J) which are non equivalent topologically. 



VI. CONCLUSIONS 



The main message of this paper was to note that there exists a systematic methods of classification and investigations 
of the FRW models with the dynamical equation of state in the quite general form p — w(a)p. It was demonstrated 
that such class of models can by reduced to a certain two dimensional dynamical system. One of the features of such 
reduction is the possibility of representation dynamics as a motion of a fictitious particle-universe in one-dimensional 
potential. The method of Hamiltonian dynamical systems can serve as a natural tool for investigation of all solutions 
for all possible initial conditions. Then one is able to see how large is the class of solutions (labelled by the inset of 
initial conditions) leading to the desired property. We consider how large is the set of initial conditions which gives 
rise to the big-rip singularity. 

In the applications of dynamics in various fields of science, the dynamics — that is dynamical system — can never by 
specified exactly. The dynamical system as a model of dynamics might be useful any way, if it can describe features 
of the phase portrait that persist when the vector field is allowed to move around. This idea emerged in the early 
history of modern dynamical system theory is called structural stability 0] . 

The idea of structural stability originated with Andronov and Pontryagin in 1937. Their work on planar systems 
was extended with the Peixoto theorem which completely characterizes the structurally stable systems on a compact, 
two dimensional manifold and establishes that they are generic. The dynamical system is said to be structurally stable 
if close dynamical systems (in some metric sense) are topologically equivalent to it. In the case of two-dimensional 
dynamical systems on the Poincare sphere there is an easy test for the structural stability of the global phase portrait 
of planar polynomial system. In particular, the phase portraits on the Poincare sphere are structurally unstable if 
there are non-hyperbolic points at infinity or if there is a trajectory connecting a saddle on the equator if the Poincare 
sphere to another saddle on S 2 p. 322). 

If we assume the existence of the moment during the evolution of the universe at which the strong energy condition 
is violated, then structural stability of the system requires the presence of saddles only (in the general case limit 
cycles and nodes are admissible). As a result we obtain three distinguished generic cases. The structural stability at 
infinity distinguishes only two phase portraits. It becomes in good agreements with phase portraits obtained from 
the reconstructed potential by using SN la 0, . 
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FIG. 1: The phase space (a, a) portrait for the FRW model with dust and dark energy. We assume the existence of the 
moment during the evolution at which p + 3p > is violated, where p and p are total energy and pressure. We have the static 
critical point marked as E at finite domain and the critical points at infinity. The system is symmetric under the reflection 
y — > —y. The trajectory of the flat model divides the phase plane on two disjoint domains occupied by closed and open 
models. Points (always static) at finite domain remains intersection of boundary of the strong energy condition and a-axis. 
The points at infinity are located on the trajectory of the flat model and circle at infinity. For large a the potential function 
looks like in the model with the cosmological constant V oc a 2 . The deSitter model deS+ is the global attractor in the future. 
The trajectories coming to this point are perpendicular to the circle at infinity S ■ The model is structurally stable because 
contains all structurally stable critical points (non- hyperbolic). 
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FIG. 2: The phase portrait for (5| (see description of Fig. in the case of V oc a m and m > 2 for large a. This is a 
representative and structurally stable model. The trajectories are tangent to the circle at infinity (da/ da = 0asa-t co). 
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FIG. 3: The phase portrait for the model under consideration (for assumption see the description of Fig. in the case of 
V(a) oc a m and m < 2 how a — » oo. In this case instead of the deS attractor the degenerated critical point a = oo, a — 
(H = 0) is present on the circle. This system is structurally unstable. 



